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Superfluid 3 He-B in the zero-temperature limit offers a unique means of studying quantum turbu¬ 
lence by the Andreev reflection of quasiparticle excitations by the vortex flow fields. We validate the 
experimental visualization of turbulence in 3 He-B by showing the relation between the vortex-line 
density and the Andreev reflectance of the vortex tangle in the first simulations of the Andreev 
reflectance by a realistic 3D vortex tangle, and comparing the results with the first experimental 
measurements able to probe quantum turbulence on length scales smaller than the inter-vortex 
separation. 

PACS numbers: 67.30.em, 47.32.C-, 47.37.+q, 67.30.lib 


Classical turbulence is well known for being simultane¬ 
ously of universal impact while analytically intractable - 
the most important unsolved problem of classical physics 
as Feynman may have expressed it. One way forward is 
to start with a simpler system. A pure superfluid in 
the T — 0 limit has zero viscosity and can be consid¬ 
ered an ideal fluid [l|. While bulk superfluid flow must 
be irrotational, it can mimic classical turbulence by sup¬ 
porting singly quantized vortices. At low temperatures, 
each vortex moves with the local superfluid velocity Q, 
comprising the combined velocity fields of all the other 
vortices [1, 3|- The system provides a concrete example 
of the vortex-filament model and the resulting complex 
flow (a vortex tangle) is quantum turbulence. 

Despite the absence of frictional dissipation, quantum 
turbulence in the T = 0 limit behaves remarkably simi¬ 
larly to classical turbulence Q exhibiting a Kolmogorov- 
like energy spectrum [E|, \&]. Turbulence in superfluid 
3 He-B at microkelvin temperatures provides several ad¬ 
vantages over other systems, most importantly that the 
vortices in this system can be visualized directly by the 
Andreev reflection of ambient thermal excitations. Such 
visualization methods have already demonstrated that a 
vortex tangle forms from the collisions of independent 
vortex rings [7], and has begun revealing statistical prop¬ 
erties of quantum turbulence d, Q. 

Here, we present the first numerical simulations of 
Andreev reflection by experimentally realistic, three- 
dimensional vortex tangles in 3 He-B, and contrast them 
with the latest experimental measurements of pure quan¬ 
tum turbulence able to probe length scales smaller than 
the average intervortex distance. This combined numer¬ 
ical/experimental approach allows us to understand the 
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connection between the vortex-tangle line density, L, (the 
total vortex-line length per unit volume) the quantity 
characterizing the intensity of turbulence, and the reflec¬ 
tion coefficient of the thermal excitations, which is used 
experimentally to visualize the turbulence. 

Andreev reflection arises in the 3 He-B Fermi superfluid 
as follows. The BCS excitation dispersion curve i£(p) 
has a minimum at the superfluid energy gap, E m i n = A, 
at the Fermi momenta pp- When an excitation transits 
from one side of a minimum to the other, its group veloc¬ 
ity reverses. For a superfluid in motion (with velocity v) 
the dispersion curve tilts by Galilean transformation to 
become E( p) + p •V 0. Thus a quasiparticle, moving 
into a region with superflow parallel to its momentum, 
experiences a potential barrier. If it has insufficient en¬ 
ergy to surmount this barrier then it must be reflected 
as a quasihole with negligible momentum transfer. Fur¬ 
thermore, most importantly, the outgoing quasihole al¬ 
most exactly retraces the path of the incoming quasiparti¬ 
cle 0. (Since any momentum transfer to the superfluid 
is minimal the excitation motion is rectilinear but the 
direction can reverse.) Similarly, quasiholes moving into 
a region of approaching flow will be Andreev-reflected 
as quasiparticles. Andreev reflection therefore offers an 
ideal passive probe for observing vortices at very low tem¬ 
peratures and can provide detailed information about the 
turbulent behavior. 

To set the scene for the simulation, we take a test vol¬ 
ume and inject a sequence of vortex rings into it. The 
rings collide, the cores intersect and recombine, gradu¬ 
ally building up an approximately homogeneous tangle. 
We then illuminate the tangle with a beam of excitations 
and calculate the reflection probability. Simulating the 
simultaneous evolution of the vortex configuration and of 
the thermal excitations is complicated and numerically 
expensive. Luckily, we can make several simplifications. 
Since the timescale of the quasiparticle motion is much 
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shorter than that of the vortex line motion {l 2 |, we first 
obtain from Eqs. © shown below the vortex configu¬ 
ration and associated flow fields, v(r, t ) at time £, and 
then analyze the propagation of excitations through this 
“frozen” flow field. Furthermore, the excitation trajecto¬ 
ries can be taken as ballistic, since the mean free paths 
at 175/iK greatly exceed any experimental dimension, 
and, finally, we assume the excitations to be point par¬ 
ticles since the coherence length, governing their spatial 
extent, ^60 nm, is tiny compared with the vortex scat¬ 
tering radius, ^20 fi m. 

The superflow field v(r, t ) and the dynamics of the 
vortex tangle are determined by the coupled equations 

V{l ' t> = -^tv^ xds ' s = v(s - t) ' (1) 

where the Biot-Savart integral extends over the entire 
vortex configuration, £, s = s(t) represents a point on the 
vortex line, and k = h/2m^ (with m 3 the mass of a bare 
3 He atom) is the superfluid 3 He quantum of circulation. 

We calculate the superfluid velocity and the time evo¬ 
lution of the vortex tangle using the vortex-filament 
method with periodic boundary conditions 0. To re¬ 
produce the experimental situation, see, e.g., Refs. 0 S, 
we take a cubic box of size D — 1 mm and numerically 
simulate the evolution of a vortex tangle generated by 
vortex loop injection for a period of 380 s. Two rings, 
radius Ri = 240 /im, are injected at opposite corners of 
the numerical domain [l4| at a frequency fi = 10 Hz. To 
ensure good isotropy, the loop injection plane is switched 
at both corners at a further slower rate f s = 3.3 Hz. The 
injected vortex loops collide and recombine, rapidly gen¬ 
erating a vortex tangle. After an initial transient of ^5 s, 
the energy content of the box reaches equilibrium. Losses 
arise from the numerical spatial resolution limit (~6 /im), 
meaning that small scale structures such as high fre¬ 
quency Kelvin waves are lost (but effectively mimicking 
the loss of kinetic energy from sound radiation at high 
frequency). The resulting tangle has an equilibrium vor¬ 
tex line density, (L) = 9.7 x 10 7 m -2 , corresponding to an 
average intervortex separation of i « (L ) -1 / 2 = 102 /im. 
The energy spectrum of this tangle is consistent with the 
k~ 5 / 3 Kolmogorov scaling for intermediate wave num¬ 
bers, fc, and with the & _1 scaling for large k , see Supple¬ 
mental Material [u| . 

To analyze the propagation of excitations, an incident 
quasiparticle flux, in (say) the x direction, is applied nor¬ 
mally to one side of the box. The quasiparticle beam is 
uniformly distributed in the (y, z) plane and covers the 
full cross section of the experimental “cell”. Ignoring an¬ 
gular factors, the incident quasiparticle flux, as a function 
of position (y, z), can be written [9] 

POO 

(™g)\v,z) = / 9{E)f{E)v g (E)dE 

J A (A ) 

= g(p F )k B T exp(-A /k B T ), 

where g(E) is the density of states, and f(JE ) is the Fermi 
distribution function, approximated at T « 0 by the 



FIG. 1. (color online). A 2D representation of the reflection 
coefficient of excitations incident on one side of the calcula¬ 
tion cell. Left: quasiparticle reflection. Right: quasihole re¬ 
flection. The vortex cores are clearly visible as the dark lines. 
The extended regions of high reflectivity (darker) and low re¬ 
flectivity (lighter) illustrate the distribution of the large-scale 
flows in the cell. 


Boltzmann distribution f(E) = exp(— E/k B T). Since 
typically quasiparticle group velocities are larger than su¬ 
perflow velocities, the quantity g(E)v g (E) in integral (|2]) 
can be replaced by g(p F ),the density of momentum 
states at the Fermi energy Q. 

In the flow field of the tangle, a quasiparticle (quasi¬ 
hole) moving against (with) a superfluid velocity v ex¬ 
periences a force dp/dt = — V(p • v), which pushes it 
towards the dispersion curve minimum where it becomes 
a quasihole (quasiparticle) with a reversed group veloc¬ 
ity. Consequently, the flux of excitations which can pass 
through a tangle is determined by the highest superfluid 
velocity, ^™ ax , encountered along the excitation’s recti¬ 
linear trajectory at constant y and 2 , and is thus given 
by: 


PC 

jYy,z = 9{PF) / 

J A 


exp {—E/k B T)dE 

! A+p F v™ ax 

= g{p F )k B T exp[—(A + p F v™ ax )/k B T \. (3) 


The fraction of quasiparticles Andreev reflected by a 
tangle along the x direction at position ( 7 /, z) is thus 


fy,z 


(nVgYy.z 

(nVg)i tZ 


1 — exp 


Pfv% 


k B T 


(4) 


The total Andreev reflection f x is the sum of the An¬ 
dreev reflections for all positions of the (y, z) plane. The 
equivalent calculation is repeated for the thermal quasi¬ 
hole flux and the results combined to give the reflection 
for a full thermal beam. 

The simulation 0 provides a large volume of informa¬ 
tion. First, we obtain the magnitude of the Andreev re¬ 
flection as a 2D contour map across the full cross section 
of the input excitation beam, see Fig. [TJ showing very 
graphically the distribution of large scale flows across 
the cell. Unfortunately, experiments do not provide us 
with similarly detailed information. Therefore, in order 
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FIG. 2. (color online). Top: the total reflection coefficient 
obtained from the numerical simulations, plotted against the 
line density of the vortex tangle. Bottom: experimental mea¬ 
surements of the fraction of Andreev scattering from vortices 
generated by a grid plotted versus grid’s velocity. Measure¬ 
ments are shown for three vibrating wires at different dis¬ 
tances from the grid as shown in the inset, see text. 


to compare theory and experiment, we concentrate in¬ 
stead on the two most significant physical outputs: the 
average Andreev reflection coefficient, (/#), and, most 
illuminating, the fluctuations of / r . 

The average calculated reflection coefficient is shown 
in Fig. [2] (top) as a function of the vortex-line density 
during the evolution of our tangle. For small line den¬ 
sities, L < 2 x 10 7 m -2 , the reflection coefficient rises 
quickly and linearly. At this stage of the tangle’s evolu¬ 
tion the injected rings are virtually noninteracting. As 
the simulation progresses, more rings enter the computa¬ 
tional domain, start to interact, collide, and form a tangle 
which absorbs all further injected rings. At higher line 
densities the rise of Andreev reflection coefficient slows, 
owing to screening effects. Here we use the term “screen¬ 
ing” to identify processes which reduce the overall reflec¬ 
tivity of the tangle for a given line density. There may be 
several mechanisms responsible and further information 
can be found in the Supplemental Material 0. 

We compare the simulation with the experiment of 


Ref. [H, which has the configuration shown in the inset of 
Fig. m The vorticity is produced by the oscillating grid 
and surrounds the vibrating wire detectors. The tangle 
flow fields Andreev reflect ambient thermal excitations 
arriving from “infinity”, shielding the wires and reduc¬ 
ing the damping. The reduction in damping on each 
wire (placed at 1.47, 2.37, and 3.49 mm from the grid) 
provides the measure of the Andreev reflection by the 
tangle, see Ref. The lower part of the figure shows 
the fractional reflection of quasiparticles incident on each 
wire. 

The numerical and experimental data plots in Fig. [2] 
have similar shapes. While the vortex line density L can¬ 
not be obtained directly from the measurements, we ex¬ 
pect the local line density of the quantum turbulence 
to increase steadily with increasing grid velocity. How¬ 
ever, the onset of turbulence is rather different in the 
simulations and in the experiment. In the simulations, 
approaching injected vortex pairs are guaranteed to col¬ 
lide and form a tangle, whereas the vibrating grid emits 
only outward-going vortex rings, with the ring flux in¬ 
creasing steadily with increasing grid velocity. At low 
grid velocities, rings propagate ballistically with few col¬ 
lisions 00. At higher velocities, ring collisions increase 
giving rise to the vortex tangle, which for ther data of 
Fig. [2] occurs when the grid velocity exceeds ^ 3mm/s. 
The data at lower velocities correspond to reflection from 
ballistic vortex rings and can be ignored for the current 
comparison. 

At higher grid velocities or tangle densities, the frac¬ 
tion of excitations Andreev reflected rises at an increas¬ 
ing rate, finally reaching a plateau. The plateau region 
is prominent in the experiments, and probably results 
from the extra quasiparticle creation produced when the 
grid reaches velocities approaching a third of the Landau 
critical velocity j22|. Compared to the simulations, the 
absolute value of the reflectivity is almost identical for the 
wire closest to the grid. This excellent agreement is per¬ 
haps fortuitous given that in the experiments quasiparti¬ 
cles travel through 1.5-2.5 mm of turbulence to reach the 
wire, compared with 1 mm in the simulation; thus, larger 
screening might be expected. A better comparison will 
require high-resolution experiments to separate the vari¬ 
ation of tangle density from the effect of the increasing 
quasiparticle numbers emitted by the grid. 

In the simulation, once the tangle has reached the sta¬ 
tistically steady state, the vortex line density and the 
Andreev reflection coefficient fluctuate around their equi¬ 
librium, time-averaged values, (L) = 9.7 x 10 7 m -2 and 
(/r) = 0.37, respectively. In order to compare the spec¬ 
tral characteristics of fluctuations of the Andreev reflec¬ 
tion 5fR(t) = /#(£) — (fit) and the vortex line density 
SL(t) = L(t) — (L) we monitor a steady state of simu¬ 
lated tangle for a period of approximately 380 s or 7500 
snapshots. Taking the Fourier transform £/#(/) of the 
time signal £/#(£), where / is frequency, we compute the 
power spectral density |£/#(/)| 2 (PSD) of the Andreev 
reflection fluctuations. Similarly we compute the PSD 
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FIG. 3. (color online). Power spectral density of the Andreev 
reflection for numerical simulations (top, blue) and experi¬ 
mental observations (middle and bottom, gray) versus fre¬ 
quency. For the details, see text. 


\SL(f)\ 2 of the vortex line density fluctuations. Fig. [3] 
shows the Andreev reflection PSD for the simulation 
(top, blue) and for the experiments (middle and bottom, 
gray). The experimental data are shown for a fully devel¬ 
oped tangle (dark gray) and ballistic vortex rings (light 
gray) relative to the grid velocities of 6.3mms _1 and 
1.9 mm s -1 , respectively. (The prominent peaks in the 
numerical data are artifacts of the discrete vortex-ring 
injection process.) 

The power spectrum of £/#(£) of the simulation and of 
the experimental data for the developed tangle (reported 
here and in Ref. [H) are in excellent agreement showing 
the same /“ 5 / 3 scaling behavior at intermediate frequen¬ 
cies. At high frequencies the experimental data develop a 
much steeper scaling (« / _3 ), not seen in the numerical 
spectrum, probably owing to the finite numerical resolu¬ 
tion. However, this frequency dependence is observed in 
the experiment where only microscopic vortex rings are 
propagating through the active region and there are no 
large-scale flows or structures. Thus, we can argue that 
the /“ 3 scaling for the vortex tangles corresponds to An¬ 
dreev reflection from superflows on length scales smaller 
than the intervortex distance. 

At a grid velocity of 6.3mms _1 , the tangle propagates 
at a mean velocity of 0.6-0.8mms _1 [9], Using Taylor’s 
frozen hypothesis, we find that the crossover between the 
two scaling laws corresponds to a length scale of ^200- 
300 fi m, in a good agreement with the intervortex dis¬ 
tance obtained from the inferred line density. 

Finally, we study the relationship between the fluctu¬ 
ations of the vortex line density, SL(t ), and the fluctua¬ 
tions of the Andreev reflection, £/#(£), by computing the 
normalised cross-correlations 


F M = + r)) 

vwrnivwim’ 


(5) 



FIG. 4. (color online). Power spectral densities of simulated 
Andreev reflection (bottom, blue) and of simulated vortex line 
density (top, red) versus frequency. Note that the accuracy 
of the initially calculated line density PSD extends to higher 
frequencies than the Andreev-reflection PSD derived from it. 
Inset: cross-correlation between the Andreev reflection coef¬ 
ficient and vortex line density. The central peak reaches a 
value of 0.9. 


where the angle brackets indicate averaging over time, 
£, in the saturated regime, and r is the time lag. The 
insert of Fig. 0] shows that the cross-correlation between 
the vortex line density and the Andreev reflection is sig¬ 
nificant, with Flr( 0) « 0.9, clearly demonstrating the 
link between them, and validating the method of visual¬ 
ization based on Andreev reflection. 

Figure [4| contrasts the spectral properties of the vor¬ 
tex line density and of the Andreev reflection from the 
simulation. At high frequencies, the vortex line spec¬ 
trum is dominated by the contribution of unpolarized, 
random vortex lines and exhibits /“ 5 / 3 scaling [l9]. In 
the intermediate frequency range, this fluctuation spec¬ 
trum shows f~ 3 scaling and is governed by the large scale 
flows indicating polarized vortex lines (polarized in the 
sense of cooperatively correlated), in agreement with re¬ 
cent numerical simulations 0- If we reasonably assume 
that this crossover should occur at around the frequency 
corresponding to the intervortex distance, £, then using 
the value of 102 fi m for i calculated above for the equilib¬ 
rium tangle, the crossover from f~ 3 to the / _5 / 3 behavior 
should occur at frequency fi ~ v/i = k/(2it£) ~ 1Hz. 
This is in very fair agreement with the frequency of the 
crossover between the two regimes of ~ 2 Hz as seen in 
Fig. m 

We conclude that the Andreev reflectance of a vor¬ 
tex tangle does indeed reveal the nature of quantum tur¬ 
bulence. The / -5 / 3 scaling of the frequency spectrum 
of the Andreev-retro-reflected signal has been observed 
earlier in the experiment of Bradley et al. ||]. Starting 
from simple physical considerations about the flow field 
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in the vicinity of a vortex filament, Bradley et al. argued 
that fluctuations of the Andreev-reflected signal can be 
interpreted as fluctuations of the vortex-line density in 
turbulent 3 He-B. The combined numerical-experimental 
results which we present here show that the fluctuations 
of the vortex line density and of the Andreev reflection 
are indeed correlated. However, their spectral densities 
behave differently. For large scale flows, the vortex line 
density scales as / -3 , while the Andreev reflection scales 
as / -5 / 3 . Interestingly, and perhaps coincidentally, the 
scaling is reversed for an unpolarized tangle. The spec¬ 


tral densities of the Andreev refection fluctuations and 
of the vortex line density scale as f~ 3 and / _5 / 3 , respec¬ 
tively. Thus it is very clear that the Andreev reflection 
technique has great potential for elucidating the proper¬ 
ties of pure quantum turbulence. 
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(1) SNAPHOT OF THE TANGLE AT THE 
MOMENT OF RING INJECTION 



FIG. 1. (Color online) A snapshot of the tangle at the mo¬ 
ment of ring injection from the opposite corners of the com¬ 
putational domain (note the vortex loops entering the com¬ 
putational box). 


(3) SCREENING MECHANISMS IN ANDREEV 
REFLECTION FROM VORTEX 
CONFIGURATIONS AND TANGLES 

The analysis of Andreev reflection from an isolated vor¬ 
tex can be found in Ref. [l|. This has been developed fur¬ 
ther in Ref. [2] to calculate the corresponding scattering 
length. 

The Andreev reflection from a single, rectilinear vortex 
line is determined by the 1/r behavior of the superflow 
field, where r is the distance from the vortex core. For di¬ 
lute vortex configurations and tangles the total Andreev 
reflection is just a sum of reflections from individual vor¬ 
tices. However, as the line density increases this is no 
longer the case due to screening effects. Here by screen¬ 
ing we mean processes which reduce the total Andreev 
reflectivity of the tangle. 


(2) ENERGY SPECTRUM OF THE VORTEX 
TANGLE 



FIG. 2. (Color online) Energy spectrum of the vortex tangle 
of the equilibrium line density ( L) — 9.7 x 10 7 m -2 . The 
slope of the spectrum is consistent with the k~ 5 ^ 3 Kolmogorov 
scaling (red dashed line) for k < kn — 27r/-£ « 6 x 10 4 m -1 , 
and with the k~ x scaling (blue dotted line) for larger k. Here 
k, m -1 is the wavenumber. 


There are two main mechanisms of screening. The 
first, which we called ‘partial screening’ in our earlier 
work M, is due to a modification of the 1/r velocity 
field of the vortex caused, in particular, by neighbour¬ 
ing vortices in its close vicinity. For example, the flow 
field at some distance from two close, (nearly) antiparal¬ 
lel vortex lines is that of a vortex dipole, v ~ 1/r 2 ; this, 
faster decay with r of the superflow velocity significantly 
reduces the total reflectivity of the vortex pair. The 1/r 
velocity field may also be modified by the effects of large 
curvature of the vortex line itself; thus, a reduction of the 
Andreev reflection from small vortex rings was reported 
in Ref. [5]. It can also be expected that the reflectivity 
of the vortex line may be reduced by the large-amplitude 
Kelvin waves. 

The other mechanism is that of ‘fractional’ screen¬ 
ing. This mechanism is of particular importance for 
large/dense tangles: vortices at the front (with respect 
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to the incident beam of excitations) of the dense or/and 
large tangle will obscure those at the rear. A simple 
model [6] of this screening mechanism was developed to 
yield a rough estimate for the vortex line density from the 
experimentally measured Andreev reflection (this mech¬ 
anism was also called ‘geometric’ screening in Ref. 0 )- 
This screening mechanism will also be important in the 


case where, in dense tangles, vortex bundles [H may form 
producing high reflectivity regions surrounded by those 
whose reflectivity is much lower. This will generate a 
significant ‘geometric’ screening since few excitations will 
reach vortices at the rear of a large vortex bundle, and 
adding a vortex to a region where there is already an in¬ 
tense reflection will not have a significant effect on the 
total reflectivity. 
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